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ABSTRACT 

• Aims. To perform a multifrequency radio study of the supernova remnant (SNR) G338.3— 0.0, in positional coincidence 

' with the TeV source HESS J1640— 465. To study the morphological and spectral properties of this remnant and its 

| surroundings searching for plausible radio counterparts to the gamma-ray emission. 

Methods. To carry out this research we observed the SNR G338.3— 0.0 using the Giant Metrewave Radio Telescope 
(GMRT) at 235, 610, and 1280 MHz. We also reprocessed archival data from the Australia Telescope Compact Array 
Oh' (ATCA) at 1290 and 2300 MHz. Also we conducted a search for radio pulsations towards a central point-like source, 

using the GMRT antennas at 610 and 1280 MHz. The molecular material in the region of the SNR was investigated 
based on observations made with the NANTEN telescope in the 12 CO (J=l-0) emission line. 

Results. The new radio observations revealed a remnant with a bilateral morphology, which at 235 MHz has the western 
wing completely attenuated because of absorption due to foreground ionized gas. The quality of the new images allows 
us to provide accurate estimates for the total radio flux density of the whole SNR at different radio frequencies. From 
the new and existing flux density estimates between 235 and 5000 MHz we derived for the whole remnant a spectral 
index a = —0.51 ±0.06 with a local free-free continuum optical depth at 235 MHz T235 = 0.9 ±0.3. No radio pulsations 
were detected towards the only radio point-like source within the HESS error circle. We derived upper limits of 2.0 and 
1.0 mjy at 610 and 1280 MHz, respectively, for the pulsed flux towards this source. No radio counterpart was found for 
the pulsar wind nebula discovered in X-rays. The inspection of the interstellar molecular gas towards G338.3— 0.0 and 
£C) . surroundings revealed that there is not any associated dense cloud that might explain a hadronic origin for the TeV 

■ detection. 
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1. Introduction ously detected by XMM-Newton and revealing a bright 

point source inside (centered at J2000 R.A.=16 h 40 m 43 s .52, 
The very high energy (VHE) 7 -ray source d ec=-46°3T35".4), labeled by the authors as Si. Based 
HESS J1640-465 was discovered by the H.E.S.S. on the mor phology and the X-ray spectrum the authors 
H , Cherenkov array during the Galactic survey 111 2004-2006 sugges t that the X-ray point source is a putative pulsar 
(|Aharonian et al. || 2006[ ). With an rms width of 2 .7 ± .5, and the diffu se componen t its associated pulsar wind neb- 
it is among the most compact sources in that sur- uk ( PWN ). |Slane et all (pOlOl) interpret the Fermi-LAT 
vey The spectrum is hard, with a photon index observations as primarily originated in such a PWN. 
T ~ 2.42, and a total integrate d flux above 200 GeV 
of ~ 2.2 x 10~ n erg cm" 2 s " 1 (lAharonian et al.1 [2006h . 

In the same area, |Slane et al.| (|2010|) reported the detec- The VHE sourcc lies in thc i ntC rior, close to the borders 

tion with the Fermi-Large Area Telescope of the source f G338.3-0.0, a shell-like radio supernova remnant (SNR). 

1FGL J 1640.8-4634 likely arising from HESS J1640-465 xhe remnailt ) about 8 ' in siz6j is located in a complex region 

Based on XMM-Newton observations, |Funk et al.| of the sky, rich in both n on-thermal and thermal emission 

(|2007[) identified a slightly extended (size ~0'.45) hard spec- ([Caswell fe Havnesll987ri . Based on HI absorption measure- 
trum X-ray emitting source (XMMU J164045.4-463131) men ts and assuming that the SNR is at the same distance 
at the centroid of the HESS source. This field has been as the neighboring HII regions, a dista nce between ~8.5 an d 
also investigated us ing the Chandra X-ray Observatory 13 kpc wa s suggested for G338 3-0.0 (jLemiere et al.ll2009h . 
flLemiere et alJl2009|) resolving the diffuse emission previ- Following iLemiere et all (|2009ft we will assume for our cal- 
culations a distance to G338.3— 0.0 of 10 kpc, leaving our 
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Up to now, the observations carried out failed to un- 
ambiguously confirm the connection between G338.3— 0.0 
and HESS J1640— 465. In this paper we present new and 
reprocessed data at 235, 610, 1280, and 2300 MHz acquired 
with the Giant Metrewave Radio Telescope (GMRT) and 
ATCA, with the aim of investigating the SNR and searching 
for a radio counterpart for the proposed pulsar and PWN 
observed in X-rays. 

2. Observations and data reduction 

2.1. GMRT observations 

We have performed low frequency observations of the 
SNR G338.3-0.0 at 235, 610, and 1280 MHz using the 
GMRT. Simultaneous dual frequency observations were 
carried out at 235 and 610 MHz during two 6 hours ses- 
sions on 26 and 27 February 2008. Data in the 1280 MHz 
band was obtained on 29 February and 1 March 2008 for a 
total 7 hours observing time on G338.3— 0.0. To make feasi- 
ble radio frequency interference (RFI) excision on the data, 
the observations at 610 and 1280 MHz were performed for 
both the upper and lower sidebands each of which with a 
total bandwidth of 16 MHz split into 128 spectral channels. 
At 235 MHz the data were collected only for the upper side- 
band using 64 channels over a total bandwidth of 8 MHz. 

At the three frequencies, the absolute amplitude scale 
was set using the calibrator 3C 286. Bandpass calibration 
was in all cases achieved by observing the source 1830—360. 
In the 235/610 MHz observations, 1625—311 was used as a 
phase calibrator, whereas 1626—298 was the secondary cal- 
ibrator for the 1280 MHz data. Automatic measurements of 
the variations in the system temperatures of the antennas at 
the GMRT were not implemented during the observations. 
Such changes in temperature may become significant when 
observing the Galactic plane region at low radio frequen- 
cies. Additional correction factors (ranged from 1.1 to 2, 
with the higher value at 235 MHz) applied to each dataset 
after imaging were estimated using interpolated values of 
sky te mperatures from the 408 MHz all-sky continuum sur- 
vey of lHaslam et al.l (119821 ) and assuming a spectral index 
/3 = —2.7 for the galactic background emission. 

The data from each day were fully reduced and im- 
aged separately to ensure that there were no day-to-day 
amplitude discrepancies. In the case of the 610 and 1280 
MHz observations, the data from each of the two avail- 
able sidebands were also reduced independently. To per- 
form this work we made use of the NRAO Astronomical 
Image Processing Software (AIPS) package following a sim- 
ilar procedure for each observed frequency. After initial cal- 
ibration the dataset from each sideband at 235, 610, and 
1280 MHz were averaged in frequency by collapsing the 
bandwidth to a number of five, eleven, and three spec- 
tral channels, respectively. The individually calibrated 235 
MHz datasets from each day were then concatenated. At 
610 and 1280 MHz, before combining data from each day, 
the data recorded for the upper and lower sidebands were 
combined into a single uv dataset using the tasks UVFLP 
and BLOAT with in AIPS and imaged as outlined below 
(|Garn et al.ll2007ft . 

In order to reduce the w-term effect imposed by the 
large field of view associated with low frequency radio ob- 
servations (the primary beam of the telescope at 610 and 
235 MHz is ^0°.8 and ~2°, respectively) we employed wide- 



field imaging as implemented in the AIPS task IMAGR. To 
deal with ionospheric based phase variations, which decor- 
relate phases across the full field of view we employed sev- 
eral rounds of phase-only self-calibration and a final phase 
and amplitude self-calibration to the visibility data at 235, 
610, and 1280 MHz. The final calibrated visibility data at 
235 MHz were combined into a single uv dataset and imaged 
using a SDI Clean algorithm. All the deconvolved images 
were finally corrected for primary beam attenuation. The 
properties of the final images are summarized in Table [TJ 

2.2. ATCA archival data 

The 1280 MHz GMRT observations are sensitive to struc- 
tures with angular scales up to 490A~ 7 arcmin. To improve 
the uv coverage we combined in the uv plane the GMRT 
data with archival ATCA data acquired in the same ra- 
dio band in November and December 2005 with the array 
operating in two different configurations (750D and 1.5C). 
The region was imaged in the mosaic mode with 5 point- 
ing centers. The ATCA data were taken simultaneously at 
1290 and 2300 MHz. For all observations, PKS 1934-638 
and PKS 1740—517 were used as primary flux density and 
secondary phase calibrator, respectively. The shortest ar- 
ray spacing of the ATCA data was 31 m corresponding to 
a maximum angular scale to which the observations were 
sensitive as large as ~ 26 and ~ 14 arcmin at 1290 and 
2300 MHz, respectively. The final angular resolution and 
sensitivity of the image at 1280 MHz formed from the com- 
bination of GMRT and ATCA data and of the image at 
2300 MHz are listed in Tabled! 

Also, to investigate in detail the radio emission around 
the X-ray pulsar candidate in a search for traces of a PWN, 
we reprocessed 4800 and 8640 MHz data corresponding to 
observations acquired with the ATCA telescope in the 750D 
configuration in 1999. These data are only useful for the 
study of structures smaller than 7 arcmin at 4800 MHz 
and 4 arcmin at 8640 MHz, which correspond to the largest 
angular scale that can be reasonably well imaged with the 
array at these high frequencies. 

3. Search for a radio pulsar 

From the archival ATCA observations carried out in the 
direction of G338. 3-0.0 at 4800 and 8640 MHz, we de- 
tected a point-like source located near the proposed X-ray 
pulsar XMMUJ164045.4-463131 dFunk et alJl2007h . unre- 
solved down to the 3" x 1" intcrferometric beam. This radio 
source, centered at R.A.=16 h 40 m 48 s , dec.=-46° 31' 58" 
(J2000), has a flux density of about 6 mJy and 4.5 mJy 
at 4800 and 8640 MHz, respectively, and a radio spectral 
index a ~ —0.5 (S cx u a ). It is worth mentioning that, 
at our current spatial resolution and sensitivity, this is the 
only radio point source present in the vicinity of the X-ray 
compact source. Though its radio spectral index is not typ- 
ical for a pulsar, separate time series observations towards 
this point source were carried out nevertheless in a search 
for radio pulses for completeness. 

The observations were conducted using the GMRT at 
610 and 1280 MHz on 24 and 25 February 2008, respec- 
tively. On both the days we used 25 antennas of GMRT 
configured in an incoherent array with a total bandwidth 
of 16 MHz. The FWHM for this configuration was 40' and 
20' for 610 and 1280 MHz respectively, centered at the point 
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Table 1. Properties of the full-resolution images of SNR G338.3— 0.0 



Frequency 


Synthesized beam 


Beam position 


Noise level 


(MHz) 


(arcsec) 


angle (deg) 


(mjy beam - ) 


235 


25.8 x 9.6 


8.0 


18 


610 


12.6 x 5.0 


-0.8 


1.6 


1280 


6.4 x 5.0 


56.9 


1.3 


2300 


4.8 x 4.3 


7.7 


0.44 



source mentioned above. This b e am al so covers the source 
Si discovered by iLemiere et al.l (|2009T > one year later. We 
also used 16 antennas located in the compact central square 
array of GMRT in a phased array mode synthesizing a more 
sensitive beam with FWHM of about 40" and 20" at 610 
and 1280 MHz respectively, centered at the point source. 
The time scries data were sampled every 256 microsec. The 
data were simultaneously acquired with incoherent array 
and phased array. This choice of observing configuration 
not only allowed us more sensitive observations of the point 
source using the phased array, but also to look for other 
related sources in the primary beam with the incoherent 
array, such as the mentioned source Si. Observations at 
610 MHz were motivated by a search for sources, which 
may be related to SNR, with a steeper spectra than the 
point source. 

The data were analyzed using the SIGPROCQ pulsar 
data analysis software on 128 processors high performance 
computer cluster currently in use at the National Centre 
for Radio Astrophysics. The data were dedispersed for 512 
trial dispersion measures (DMs) with DMs ranging from 
to about 465 pc cm -3 at 610 MHz. DMs were spaced more 
coarsely at 1280 MHz, due to the smaller dispersion smear- 
ing across each individual frequency channel. At this fre- 
quency the dedispersion search was made in 160 steps using 
DMs between to 660 pc cm -3 . We searched for periodici- 
ties in the dedispersed time series using a standard harmon- 
ics search and used a threshold signal-to- noise (S/N) of 8. 
Interference periodicities corresponding to power line RFI 
as well as other known RFI were eliminated. After a careful 
examination of the profiles for the candidates produced by 
search processing, we did not find any characteristic peri- 
odic signal, as expected from a pulsar. 

A calibrator source, 1830—360, was also observed using 
both the incoherent and phased array. This source was used 
to estimated the rms noise in the time series. Upper limit 
on flux density of any pulsar in our field of view with S/N 
of 8 times rms were estimated as 2 and 1 mjy for 610 and 
1280 MHz, respectively. Therefore, though our observations 
did not reveal pulsations from an intense source, they do 
not rule out a pulsar fainter than the limits quoted above. 
In addition, the pulsed radiation may be beamed away from 
our line of sight. 

With a spectral index of —0.5, much less steeper 
than the average pulsar spectral index of about —1.8 
(jMaron et al.ll2000t) . our limits on pulsed radiation imply 
that the radio point source, shown in Fig. [1] must be un- 
related. Source Si may be a radio pulsar, but is either 
fainter than our limits or beamed away. Lastly, the high 
sky background and large interstellar scattering at a dis- 
tance of 10 kpc could also lead to non-detection of pulsed 



1 More information available at |http://sigproc.sourceforge.net 



radiation, particularly for periods below 100 ms, due to 
scattcr-broadening of the radio pulse. 

4. New radio images of G338.3— 0.0 

In Fig. [T] we present close-up views of a portion of the 
field mapped with GMRT at the low radio frequencies plus 
the ATCA reprocessed data, showing the detailed mor- 
phology of G338.3-0.0 at 235, 610, 1280, and 2300 MHz. 
These are the first high-resolution images of the source 
providing unprecedented sensitivity to the radio emission. 
The fact that the smallest baseline in the GMRT Central 
Square for which the data at 235 and 610 MHz were reli- 
able correspond to ~ 86A and ~ 230A respectively, equiva- 
lent to structures with angular scales up to 40 and 15 ar- 
cmin, ensures that the contributions from the largest spa- 
tial scales have been adequately recovered for this SNR 
about 8' in size. Particularly, the GMRT 610 MHz im- 
age (Fig. [T]d) improves by a factor 6 the angular resolu- 
tion and 7 times the sensitivity in comparison with the 
only previously known im age of this remnant at 843 MHz 
(IWhiteoak fc Greenl Il996h . In the case of the data at 
1280 MHz, since G338.3-0.0 is almost at the limit of 
the largest imaged structures we incorporated the shortest 
spacings information as taken from archival ATCA observa- 
tions to avoid missing large scale flux density in the image 
at this frequency, as described in Sect. 2.2. 

Except for the radio image at 235 MHz, where most of 
the western half is completely attenuated, the total inten- 
sity images at 610, 1280, and 2300 MHz are quite similar 
in appearance, revealing a supernova remnant with a clear 
bilateral morphology. Considerable small scale features, in- 
cluding bright knots, are evident in the new images along 
the entire shell, with their respective brightness varying 
across the spectrum. Diffuse emission is also observed in- 
side the radio shell, particularly close to the eastern radio 
shell. 

The large sky field mapped at 610 MHz allows us 
to get a picture of the different sources emitting around 
G338.3— 0.0. In Fig. [2] we present the radio emission at 
610 MHz in an area of about 0.5 square degrees. Several 
HII regions are evident in the field. The position and ex- 
tension of the cataloged thermal regions (from Caswell & 
Hayncs 1987 and Paladini et al. 2003) arc indicated in Fig. [5] 
by white dashed circles. From this figure it is obvious that 
the remnant lies in a complex region in our Galaxy with 
an intricate network of non-thermal and thermal structure, 
with ample regions of the SNR shell covered by the HII 
regions G338.45+0.06 and G338.4+00.0. 

To accurately establish the flux density and hence the 
spectral properties of the SNR emission it is very important 
to disentangle thermal from non-thermal contributions. To 
carry out this task we made use of the Spitzer GLIMPSE 
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Fig. 1. New radio continuum images from the SNR G338.3— 0.0. The contours and the linear intensity scales have been 
chosen to emphasize the relevant features at each frequency. The wedges display in mJy beam -1 the values of the radio 
emission at each frequency. All images were aligned to the same pixel position. The coordinates displayed in Fig. lb 
are the same for the other three panels, a) GMRT image at 235 MHz. The emission towards the western hemisphere 
is significantly suppressed at this frequency. The contours are traced at 36, 50, 60, 74, 90, and 120 mJy beam -1 , b) 
GMRT image at 610 MHz. The intensity contours correspond to 3.5, 5, 6, 7, 8.5, and 10 mJy beam -1 , c) Radio emission 
at 1280 MHz from GMRT and ATCA data. The contours are traced at 4, 5, 6, 8, and 10 mJy beam -1 , d) 2300 MHz 
image obtained from reprocessed archival ATCA data. The radio contours are traced at 2, 3, 4, and 6 mJy beam -1 . All 
the displayed images include primary beam correction. The synthesized beams, position angle and noise levels are listed 
in Table [TJ The white dashed circle included in Fig. 1 (b) shows the position and size of the VHE gam ma-ray source 
HESS J1640— 465. The plus sign in each panel marks the position of the X-ray source Si identified by iLemiere et al.l 
(|2009h . (A color version of this figure is available in the on line journal). 



and MIPSGAL infrared data. In Fig. [3] we present a di- 
rect comparison of the IR emission with the new GMRT 
610 MHz image in the spatial region around G338.3— 0.0. 
The emission in the 8 and 24 ^m infrared bands are ex- 
cellent tracers of warm dust associated with star forming 
regions, usually related to thermal radio emission. 

Figure [3] clearly demonstrates that the entire complex 
of HII regions located to the north of G338.3— 0.0 with an 
extent of about 12' in the sky, is prominent in radio and 
in the two IR bands. In addition to the already noticed 
thermal emission, two infrared condensations particularly 
bright in 24 fj,m, arc observed overlapping the SNR. One 
located on the southern extreme of the western shell and 
the other near the center of the eastern half of the remnant. 
These sources with both strong infrared and radio emis- 



sion are probably normal or ultracompact HII regions over- 
lapping the non-thermal SNR emission. Other noticeable 
structure is a bright ring of infrared 8 /im and 24 fj,m emis- 
sion located on the northwestern portion of the SNR, which 
appears to outline precisely an enhancement in the radio 
band. Suc h a feature was fir s t iden tified at infrared wave- 
lengths by IChurchwell et all (|2006l ) using GLIMPSE data 
who proposed that it is the projection of a fully enclosed 
three-dimensional bubble (designated as S33 in Churchwcll 
et al.'s catalog) with an angular diameter of about 0.9 ar- 
cmin. The composite image also reveals the presence of faint 
infrared emission at 8 /im along the northwestern boundary, 
within the error circle of the VHE source HESS J1640— 465. 
We also investigated the near infrared emission using the 
3.6 /im, 4.5 /im, and 5.8 /im Infrared Array Camera (IRAC) 
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16h40m00.0s 30 0s 



Fig. 2. Radio continuum image at 610 MHz of a large area 
around the SNR G338. 3—0.0, showing several nearby ther- 
mal sources encircled by dashed white lines to facilitate 
their location. The intensity scale is based on a linear re- 
lation from 2.4 to 8 mJy beam -1 . (A color version of this 
figure is available in the on line journal). 



data from the GLIMPSE Science Program but no emission 
above noise is detected in this region. 




With this information we attempt to provide ac- 
curate estimates for the flux densities over the whole 
SNR G338. 3-0.0 at 235, 610, 1280, and 2300 MHz. In ad- 
dition we used the image of this SNR acquired at 843 MHz 
with the MOST telescope and the image at 5000 MHz taken 
from Parkes 64-m telescope archival data to extract another 
values for the flux density. Our estimate of the flux density 
at 843 MHz is in good agreement, wit h in un certainties, with 
that reported bv lWhiteoak fc Greenl (|1996l ). 

These flux density estimates are corrected by a back- 
ground level that was determined at each frequency by 
tracing one-dimensional plots of the intensity as a func- 
tion of position in selected slices around the shell of the 
SNR. We note a dependence on the direction of measure- 
ment of the background emission specially for the images 
at 1280 and 2300 MHz, for which variations of a few tens 
of mJy beam -1 were calculated between the east and west 
directions of the radio shell. The values estimated in this 
way were then subtracted from the integrated flux density 
at each frequency. An extra contribution from the thermal 
emission was also subtracted from the integrated flux den- 
sity measurements. This contribution was estimated based 
on the comparison of our radio image with the infrared data 
by integrating the infrared flux in the bright regions over- 
lapping the SNR and determining their respective fluxes at 
the observed radio frequencies by assuming a spectral index 
of about a ~ —0.1 for the thermal emission. Table [21 gives 
the integrated flux density measurements for G338.3— 0.0 
from 235 to 5000 MHz. In the second column we list the re- 
sults obtained by directly measuring the flux density from 
the image at each frequency. The listed flux densities rep- 
resent the average from different measurements carried out 
using slightly different contours to enclose the SNR emit- 
ting region. In columns 3 and 4 we summarize the back- 
ground and thermal emission contribution, while the final 
flux density values after subtracting these contributions are 
listed in the fifth column of Tablc[2] The results in this table 
have been p l aced o n the same absolute flux density scale of 
iBaars et all (|1977t) . 

Figure 2] shows the radio continuum spectrum of 
G338.3— 0.0 from the flux densities summarized in Table [5J 
There is a clear sign of a low frequency turnover in this spec- 
trum likely to be caused by the emission associated with the 
intervening thermal gas. We first used a single power law 
spectrum excluding the flux density value at 235 MHz to de- 
termine the integrated spectral index of the remnant (rep- 
resented as a dashed line in Fig. @|. A weighted fit yields a 
spectral index a = —0.47 ±0.05 (S oc v a \ typical for shell- 
type SNRs. In addition, to determine the properties of the 
gas whose emission is responsible for the observed turnover, 
we have used all the integrated flux measurements, to make 
a weighted fit to the equation (indicated by the solid line 
in Fig. |U) 



Fig. 3. A comparison of the radio continuum emission at 
610 MHz (in blue) from the region around SNR G338.3-0.0 
and the infrared emission at 8 fiia (in green) and 24 /im (in 
red) taken from the Spitzer Space Telescope GLIMPSE and 
MIPSGAL surveys, respectively. The white circle marks 
the position and extension of HESS J1640— 465, while 
the cross indicates the location of the X-rays source Si 
(jLemiere et al.ll2009t ). (A color version of this figure is avail- 
able in the on line journal) . 



S » = Si0S (i08lHfe) QeXP - T408 ( 



,408 MHz 



-2.1 



(1) 

where 5 y (Jy) represents the flux density at the frequency 
f(MHz). S408 and T408 are the flux density and the aver- 
age optical depth at the fiducial frequency of 408 MHz. 
This equation assumes a constant spectral index a for the 
synchrotron emission over the entire radio band, a l lowin g 
for a turnover at the lowest frequencies (|Kassimlll989t) . 
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Table 2. Flux density estimates for the SNR G338.3-0.0 



Frequency 
(MHz) 


Directly measured 
flux density (Jy) 


Background 
emission (Jy) 


Thermal 
emission (Jy) 


Final 
flux density (Jy) 


235 


7.0 ± 1.2 


0.5 ±0.2 


0.3 ±0.2 


6.2 ± 1.2 


610 


8.80 ± 0.94 


0.20 ±0.05 


0.45 ±0.10 


8.15 ±0.95 


843 


8.6 ±0.7 


0.3 ±0.1 


1.2 ±0.4 


7.1 ±0.8 


1280 


8.5 ±0.7 


0.9 ±0.3 


1.3 ±0.2 


6.3 ±0.8 


2300 


7.1 ±0.7 


1.1 ±0.4 


1.2 ±0.1 


4.8 ±0.8 


5000 


4.7 ±0.2 


0.3 ±0.1 


1.3 ±0.1 


3.1 ±0.2 




to the western and northwestern part of the SNR where T235 
varies between 1 and 2. Based on this result together with 
our spectral analysis, we can derive the physical properties 
of the thermal plasma responsible for the observed atten- 
uation in the western part of the SNR. We conclude that 
if the electron temperature varies in the typical range for 
HII regions (5000 to 10000 K) and if the absorbing gas has 
a deepness along the line of sight comparable to its trans- 
verse size, at the assumed distance of 10 kpc, the emission 
measure varies between 0.1 and 0.3xl0 6 pc cm -6 in the 
region where the absorption is stronger, and the electron 
density is in the range 100 to 165 cm~ 3 . 



0.2 r i 

2.0 2.5 3.0 3.5 4.0 

Log v (MHz) 

Fig. 4. Integrated radio continuum spectra of the 
SNR G338.3— 0.0 obtained from the flux density values 
listed in Tabled (column 5). The dashed line is the best 
fit to the flux density measurements at frequencies higher 
than 235 MHz, which gives an integrated spectral index 
a = -0.47 ± 0.05. The solid line is the fit to the data 
using Eq.([T]), which takes into account a low frequency 
turnover due to absorption by foreground thermal gas. 
This method yields a global spectral index for G338.3— 0.0 
a = -0.51 ±0.06. 

The best fit for the parameters are 5*408 = 11.2 ± 1.3 Jy 
and a = —0.51 ± 0.06, whereas the best fit value for the 
free-free continuum optical depth is T408 = 0.28 ± 0.08. 
The average optical depth at 235 MHz estimated from 
t v = t 408 [^/408]~ 21 is 0.9 ± 0.3, which should be consid- 
ered as an upper limit due to the fact that the low frequency 
turnover observed in the spectrum of Fig. U occurs near the 
lowest o bserved frequ ency, resulting in a possible overesti- 
mation |Kassim 1 9891) . Accurate flux density measurements 
of G338.3— 0.0 at very low frequencies (< 100 MHz) are 
highly desirable to better establish the frequency depen- 
dence of the spectrum in this region. The derived spectral 
properties of G338.3— 0.0 are not surprising given that, as 
we have shown in Fig. [5J the remnant lies in a site that 
hosts numerous HII regions. Thus, contamination from ion- 
ized gas from these HII regions (or their envelopes) along 
the line of sight constitutes a potential candidate to provide 
extrinsic free-free absorption towards the SNR. 

By applying appropriate uv-tapcring, in order to match 
the range of spatial scales measured at 235 and 610 MHz, 
we studied the dependence on the position across the face 
of the remnant of the optical depth at the lowest frequency. 
As expected the higher values of this magnitude correspond 



5. Search for a PWN in G338.3-0.0 

Based on the new radio images we have searched for a 
counterpart to the observed X-ray emissi on interior to the 
SNR radio shell, pro posed to be a PWN (|Funk et al.ll2007l: 
iLemiere et al.ll2009|) . However, down to the noise level of 
our data at each frequency they do not reveal any trace 
of radio emission associated with the X-ray PWN. On the 
contrary, instead of the expected nebular enhancement, a 
depression in the emission is observed at 610 MHz in the 
area covered by the X-ray nebula. In spite of the complex 
nature of the emission in this region, there is no reason to 
suspect a deficit in the recovered flux density, since obser- 
vations at this frequency were made with very good spatial 
coverage and no overlap of thermal emission is evident in 
coincidence with the X-ray nebula. We can therefore only 
set limits on the radio flux density assuming for our cal- 
culation that particles originate close to the X-ray point 
source and then disperse in a nebula whose extent in ra- 
dio is, at least, comparable to that of the X-ray nebula in 
G338.3— 0.0. The best constraint is ~ 6 mJy obtained from 
the image at 610 MHz as it represents the lowest between 
all the flux density measurements. 

Assuming a PWN radio spectral index of a = —0.3 
(jGaensler et al.l l2000t) , the mentioned upper limit corre- 
sponds to a broadband radio luminosity integrated between 
10 8 and 10 12 Hz of L n ~ 6 X 10 32 d 2 erg s~\ where the 
upper cutoff frequency considered in this calculation is ob- 
tained from the spectral break betw een radio and X- rays 
from the XMM-Newton observations (jFunk et al.ll2007l ). In 
this case, the derived radio luminosity for this object is dis- 
tinctly lower than Lr ~ 10 34 erg s" 1 , observed in other 
radio PWNe (jGaensler fc Sland 120061 ) . On the contrary, a 
higher luminosity of Lr ~ 9.0 x 10 34 d 2 erg s _1 is ob- 
tained if, the upper cuto ff is derived from th e fit between 
radio and Chandra data (jLemiere et al. 2009J, i.e. a lower 
and upper cutoff frequency of 10 8 and 10 16 Hz are used in 
the calculation. This latter result, however, is at odds with 
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the non-detection of a radio nebula in the region. The dra- 
matic difference between the calculated luminosities comes 
from the uncertainty in the frequency at which the spectral 
break between radio and X-rays takes place, which is due 
to the quite different X-ray slopes quoted in th e literature. 
In ord er to reproduce the broad-band emission iFunk et al.l 
(|2007t ) proposed a time-dependent rate injection of rela- 
tivistic electrons in which the X-ray emitting region is orig- 
inated by youngest particles (t < 2000 yr), while an older 
population of electrons injected early in the formation of the 
PWN dominates the VHE 7-ray emission. In this scenario 
if the population of youngest electrons also dominates the 
radio emission, the values of the radio flux densities would 
fall well below the sensitivity limit of current observations, 
compatible with our non-detection of the radio PWN. 

6. The environmental molecular gas 

In a search for an alternative explanation for the observed 
7-ray emission, we analyzed the 12 CO (J=l— 0) molecular 
line emission around G338.3— 0.0 based on data taken from 
the CO survey carried out with the 4-m NANTEN telescope 
operated in Las Campanas Observatory (Chile). These ob- 
servations have a HPBW of 2.7 arcmin, a grid spacing of 
4 arcmin, and a velocity resolution of 1 km s _1 . 

After carefully inspecting the complete 12 CO data cube 
we did not find any molecular cloud that, based on morpho- 
logical correspondence, might be associated with either the 
source HESS J1640-465 or the SNR radio shell, at least 
down to the sensitivity and resolution of the CO survey 
used. Deeper molecular observations would be very useful 
to investigate the likelihood of a hadronic origin for the 
observed 7-ray emission. 

7. Summary 

This paper presents the highest angular resolution and 
sensitivity multifrequency dataset of the radio continuum 
emission towards the SNR G338.3— 0.0. A mixture of syn- 
chrotron radiation from the remnant and overlapping unre- 
lated thermal gas it is noticeable in the mapped region. The 
quality of the new images provides for the first time proper 
estimates for the total radio flux density of the whole SNR 
at four radio frequencies. Based on these images we stud- 
ied the global radio continuum spectrum of the remnant 
and found a clear indication of a low frequency turnover. 
We conclude that a foreground sheet of ionized gas is ab- 
sorbing the radio synchrotron emission coming from almost 
half of the radio SNR shell. We interpret that this absorp- 
tion is associated with the complex of HII regions in close 
proximity to the remnant. 

We also identified a previously uncataloged point- 
like radio source near the centroid of the TeV source 
HESS J 1640— 465, close to t he compact object detected in 
X-rays (jLemiere et al.l [20091 ) and considered by these au- 
thors as a putative pulsar which would be powering a PWN 
detected in X-rays. No pulsed emission associated with the 
radio point source, or within the field of view, was detected 
up to a continuum flux density of 2.0 and 1.0 mJy at 610 
and 1280 MHz, respectively. Besides, no nebular radio emis- 
sion was detected in correspondence with the X-ray PWN 
in spite of the good quality of our radio images down to 
low surface brightness limits. In this context, we only set 



an upper limit on the radio emission at 610 MHz over a 
region with an area comparable to that of the PWN seen 
in X-rays. 

Finally, we analyzed the environmental conditions for a 
hadronic origin of the 7-ray emission. We concluded that 
no bright molecular matter is apparent overlapping cither 
the source HESS J1640-465 or the SNR radio shell. 

In summary, these radio observations fail to provide fur- 
ther support to the association of the SNR and/or its PWN 
with the emission observed in the 7-rays range. Physical 
mechanism producing only high energy photons should be 
proposed to reconcile the observations across the electro- 
magnetic spectrum. 
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